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O '• 

CNJ ■ Dyson-Schwinger equations furnish a Poincare covariant approach to hadron physics. 

They reveal that dynamical chiral symmetry breaking is tied to the long-range behaviour 
of the strong interaction and make predictions corroborated by modern lattice-QCD simu- 
lations. A hallmark in the contemporary use of DSEs is the existence of a nonperturbative, 
symmetry preserving truncation that enables the proof of exact results. The systematic 
error associated with the truncation's leading term has been quantified and this underpins 
an efficacious one-parameter model that is being employed to study meson excited states. 

in : 

1. DRESSED-QUARK PROPAGATOR 

The Dyson-Schwinger equations (DSEs) are a nonperturbative approach to studying 
continuum QCD [ 12] • At the simplest level they provide a generating tool for per- 
turbation theory and, because QCD is asymptotically free, this feature has the potential 
to materially reduce any model dependence in DSE applications. This comment is ex- 
plained by noting that the long-range behaviour of the interaction between light-quarks is 
hitherto unknown and must therefore be modelled. Hence, the fact that a weak coupling 
expansion of the DSEs reproduces every diagram in perturbation theory entails that in 
>• . this approach the model dependence is restricted to the infrared domain k 2 < 1 GeV 2 . 

The DSEs are a nonperturbative tool and can therefore be used directly to study: 
hadrons as bound states; the importance of dynamical chiral symmetry breaking (DCSB); 
and the confinement of quarks and gluons. These phenomena are all linked to the infrared 
behaviour of the interaction. Hence the DSEs provide a bridge between this behaviour 
and observables, and thereby a means by which it can be charted. 

The best known DSE is the simplest: the gap equation, which describes how a fermion's 
propagator is modified by interactions with the medium being traversed. In QCD that 
equation is the DSE for the dressed-quark propagator: 

Sip)' 1 = Z 2 (i 1 -p + m barc ) + Z x J* g 2 D^(p - q )^S(q)T a v {q; p) , (1) 

wherein: D^ y {k) is the dressed-gluon propagator; T"(g;p) is the dressed- quark-gluon ver- 
tex; m bare is the A-dependent current-quark bare mass; and / A := J A <i 4 g/(27r) 4 repre- 
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sents a translationally-invariant regularisation of the integral, with A the regularisation 
mass-scale. In addition, Z 1)2 (C 2 , A 2 ) are the quark-gluon- vertex and quark wave function 
renormalisation constants, which depend on A and the renormalisation point, (, as does 
the mass renormalisation constant Z m (( 2 ,A 2 ) = Z 4 (( 2 , A 2 )/Z 2 (( 2 , A 2 ). The solution has 
the form 

S-\ P ) = %1 . P A( P 2 , C 2 ) + B(p 2 , C 2 ) = 1 [i!-p + M(p 2 )} , (2) 

where M(( 2 ) = m(Q '■= ^bare(A) Z~ 1 (^ 2 , A 2 ) is the running quark mass. 

The gap equation illustrates the features and flaws of each DSE. It is a nonlinear 
integral equation for S(p) and hence can yield much- needed nonperturbative information. 
However, the kernel involves the two-point function D^k) and the three-point function 

(<?;£>)• The gap equation is therefore coupled to the DSEs these functions satisfy. 
Those equations in turn involve higher n-point functions and hence the DSEs are a tower 
of coupled integral equations with a tractable problem obtained only once a truncation 
scheme is specified. It is unsurprising that the best known truncation scheme is the weak 
coupling expansion, which reproduces every diagram in perturbation theory. This scheme 
is systematic and valuable in the analysis of large momentum transfer phenomena because 
QCD is asymptotically free but it precludes any possibility of obtaining nonperturbative 
information, which we identified as a key aspect of the DSEs. 

In spite of the problem with truncation, gap equations have long been used efficaciously 
in obtaining nonperturbative information about many-body systems. The positive out- 
comes have been achieved through the simple expedient of employing a rudimentary 
truncation and comparing the results with observations. Naturally, agreement under 
these circumstances is not compelling evidence that the contributions omitted are small. 
However, it does justify further study, and an accumulation of good results is grounds for 
a concerted attempt to substantiate a reinterpretation of the truncation as the first term 
in a systematic and reliable approximation. That has recently been achieved 

The gap equation's solution has long been of interest in grappling with DCSB and 
typical results are depicted in Fig.^ The infrared suppression of Z(p 2 ) and enhancement 
of M(p 2 ) evident in this figure are longstanding predictions of DSE studies [H], and in this 
behaviour one perceives directly the evolution of a current-quark into a constituent-quark. 
A critical feature is that, so long as the kernel of the gap equation has sufficient integrated 
strength on the infrared domain, a nonzero running quark mass is obtained even in the 
chiral limit. This effect is DCSB. It is impossible at any finite order of perturbation 
theory and also apparent in Fig. ^ The precise connection between the dressed-quark 
propagator and vacuum quark condensate is discussed in Ref. and a connection with 
quark confinement is canvassed in Ref. 

The dressed-quark propagator can be calculated in lattice-regularised QCD. Results are 
available in the quenched truncation, and depicted in Fig.[T]are those of Ref. [Ej obtained 
with the current-quark masses (£ = 19 GeV) 

(3) 



a ^lattice 


0.018 


0.036 


0.072 


m(C)(GeV) 


0.030 


0.055 


0.110 



The precise agreement with DSE results is not accidental. (NB. The chiral limit discrep- 
ancy points to the inadequacy of a linear extrapolation of lattice data to this limit 



3 




0.0 1.0 2.0 3.0 4.0 0.0 1.0 2.0 3.0 4.0 

p (GeV) p (GeV) 

Figure 1. Left Panel - Dashed-curve: gap equation's solution in the chiral limit; solid 
curves: solutions for M(p 2 ) obtained using m(() in Eq. (jHJ). (From Ref. Data, upper 
three sets: lattice results for M{p 2 ) in GeV at am values in Eq. (jSJ); lower points (boxes): 
linear extrapolation of these results [|7] to am = 0. Right Panel - Dashed curve, Z(p 2 ), 
and solid curve, M(p 2 ) calculated from the gap equation with m(() = 55MeV [E]. Data, 
quenched lattice-QCD results for M(p 2 ) and Z(p 2 ) obtained with am = 0.036 [|7j. 



The essential agreement between lattice results and DSE predictions was highlighted in 
Refs. [IEJ E] but Ref. [Ej pursued a different goal. Only recently has reliable information 
about the gap equation's kernel at infrared momenta begun to emerge, in the continuum 
[ QUI and on the lattice [ • Reference [ HJ| therefore employed an Ansatz for the in- 
frared behaviour of the gap equation's kernel in order to demonstrate that it is possible 
to correlate lattice results for the gluon and quark Schwinger functions via QCD's gap 
equation. This required the gap equation's kernel to exhibit infrared enhancement over 
and above that observed in the gluon propagator alone, which could be attributed to an 
amplification of the dressed-quark-gluon vertex whose magnitude is consistent with that 
observed in quenched lattice estimates of this three-point function [I12j. 

2. HADRONS 

It is evident that reliable knowledge of QCD's two-point functions (the propagators for 
QCD's elementary excitations) is available. Direct comparison with experiment requires 
an equally good understanding of bound states. Progress here has required the evolution 
of an understanding of the intimate connection between symmetries and DSE truncation 
schemes. This is well illustrated by considering the pion, whose properties are profoundly 
connected with DCSB. Indeed, the correct understanding of pion observables requires an 
approach to contain a well-defined and valid chiral limit. 

Chiral symmetry and its breaking are expressed through the axial-vector Ward-Taka- 
hashi identity 

PpTi^k- P) = ^(A^y + i l5 T -^S-\k-) - iM(C) n(k; P) ~ Ti(k; P) iM(C), (4) 

where k± = k ± P/2, and {r- 7 , j = 1,2,3} are the Pauli matrices as herein we focus 
on S77(2)-flavour. This identity connects the axial-vector vertex: F 3 5 Jk;P), P is the 
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Figure 2. Planar dressed-quark-gluon vertex obtained by neglecting contributions associ- 
ated with explicit gluon self-interactions. Solid circles indicate fully dressed propagators. 
The quark-gluon vertices are not dressed. 



total momentum, with the dressed quark propagator: S = di&g[S u , Sd), the pseudoscalar 
vertex: T J 5 (k;P), and the current-quark mass matrix: -M(C) = diag[m u (£), md(()]- The 
propagator satisfies the gap equation but the vertices are determined by inhomogeneous 
Bethe-Salpeter equations; e.g., 



1 5^ 



III 









tu h 







[S(q + )T^(q; P)S(q.)} sr K£(q, fc; P) 



(5) 



wherein K(q, k\ P) is the dressed-quark-antiquark scattering kernel. The importance of 
DCSB entails that any truncation useful in understanding low energy phenomena must 
be nonperturbative and preserve Eq. (0} without fine tuning. This constraint requires an 
intimate connection between K(q, k; P) and the gap equation's kernel. 

One systematic and nonperturbative truncation scheme has been identified that expli- 
cates this connection, and hence preserves QCD's global symmetries [0]. It is a dressed- 
loop expansion of the dressed-quark-gluon vertices that appear in the half-amputated 
dressed-quark-antiquark scattering matrix: S 2 K. In this scheme, as in perturbation the- 
ory, it is impossible, in general, to obtain complete closed-form expressions for the kernels 
of the gap and Bethe-Salpeter equations. However, for the planar dressed-quark-gluon 
vertex depicted in Fig. El closed forms can be obtained and a number of significant features 
illustrated when one uses the following model for the dressed-gluon line [ FEJj 



g'D^ik) 



k^ky 



J /1V 



(27Tfg 2 5\k) 



(6) 



where Q sets the model's mass-scale. This form has many positive attributes in common 
with the class of renormalisation-group-improved models and furthermore its particular 
momentum-dependence works to advantage in reducing integral equations to character- 
preserving algebraic equations. 

It is a general result that with any vertex whose diagrammatic content is explicitly 
known; e.g., Fig.El it is always possible to construct a unique Bethe-Salpeter kernel 
which ensures the Ward-Takahashi identities are automatically fulfilled. That kernel is 
necessarily non-planar. This becomes transparent with the model in Eq. (jHJ), using which 
the gap equation obtained with the vertex depicted in Fig.|21 reduces to an algebraic 
equation, irrespective of the number of dressed-gluon rungs that are retained, and the 
same is true of the Bethe-Salpeter equations in every channel: pseudoscalar, vector, etc. 
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Tabic 1 



7T and p meson masses obtained with Q = 0.48 GeV. (Dimensioned quantities in GeV.) n 
is the number of dressed-gluon rungs retained in the planar vertex, see Fig. |2l and hence 
the order of the vertex-consistent Bethe-Salpeter kernel. 
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M n=oo 


7r, m = 














7r, m = 0.011 


0.152 


0.152 


0.152 


0.152 


p, m = 


0.678 


0.745 


0.754 


0.754 


p,m = 0.011 


0.695 


0.762 


0.770 


0.770 



Results for the it and p are illustrated in Table H It is evident that, irrespective of 
the order of the truncation; i.e., the number of dressed gluon rungs in the quark-gluon 
vertex, the pion is massless in the chiral limit. (NB. This pion is composed of heavy 
dressed- quarks, as is evident in the calculated scale of the dynamically generated dressed- 
quark mass function: M(0) ~ 0.5 GeV; viz., dynamically generated constituent-quarks 
compose the pion.) The masslessness of the 7r is a model-independent consequence of 
the consistency between the Bethe-Salpeter kernel and the kernel in the gap equation. 
Furthermore, the bulk of the p-7r mass splitting is present in the chiral limit and with the 
simplest (n = 0; i.e., rainbow-ladder) kernel, which makes plain that this mass difference 
is driven by the DCSB mechanism: it is not the result of a finely adjusted hyperfine 
interaction. Finally, the quantitative effect of improving on the rainbow-ladder truncation; 
viz., including more dressed-gluon rungs in the gap equation's kernel and consistently 
improving the kernel in the Bethe-Salpeter kernel, is a 10% correction to the vector meson 
mass. Simply including the first correction (n = 1; i.e., retaining the first two diagrams 
in Fig. |2I) yields a vector meson mass that differs from the fully resummed result by ~ 1%. 
The rainbow-ladder truncation is clearly accurate in these channels. 



3. AB INITIO CALCULATIONS 



The leading-order term in the truncation scheme described above is the renormalisation- 
group-improved rainbow- ladder truncation, which as we have just seen is reliable for 
flavour nonsinglet pseudoscalar and vector mesons. This truncation has been refined and 
exploited by Maris and Tandy in a series of articles 1 via a one-parameter model for the 
behaviour of the dressed-quark-quark interaction: 



Q 2 



- 6 DQ 2 e 



+ 



iln 



T + (l + Q 2 /A 2 



QCD 



HQ 2 ) 



(?) 



wherein F{Q 2 ) = [1 - exp(-Q 2 /[4m 2 ])]/Q 2 , m t = 0.5 GeV; r = e 2 - 1; 7m = 12/25; 
and Aqcd = 0.234 GeV. This simple form expresses the interaction as a sum: the second 
term ensures that perturbative behaviour is preserved at short-range; and the first makes 
provision for enhancement at long-range. The true parameters in Eq. (JJJ) are D and u, 



1 An exemplary success was their prediction [ I14| of the electromagnetic pion form factor [I15j. In addition, 
see Ref. [El an d references therein, and the review of their substantial contributions in Ref. 
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which together determine the integrated infrared strength of the rainbow-ladder kernel. 
However, they are not independent: in fitting to a selection of observables, a change 
in one is compensated by altering the other; e.g., on the domain uj e [0.3, 0.5] GeV, the 
fitted observables are approximately constant along the trajectory [[8] u) D = (0.72 GeV) 3 . 
Hence Eq. (jTJ) is a one-parameter model. This correlation: a reduction in D compensating 
an increase in u, ensures a fixed value of the integrated infrared strength. It also defines a 
single dressed-glue mass-scale, m g = 720 MeV, which characterises infrared gluodynamics. 

A new application is the calculation of properties of meson excited states, which is a 
first step in seeking exotics and hybrids. As an example, it is a general feature of QCD 
that the axial- vector vertex, obtained via Eq. (|5|). exhibits a pole whenever P 2 = —m 2 n , 
where m Wn is the mass of any pseudoscalar u-d meson; viz., 

r ^ fc; p )U m?r = regular terms + P 2+^ r i^; p ) > ( 8 ) 

with r{ n (/c; P) being the 0~ + bound state's Bethe-Salpeter amplitude: 

Ti n (k;P) = r j l5 [iE^(k;P) + 1 -PF 7rn (k;P) 

+ 1 -kk-PG 7rn {k;P) + a tlu k fl P u H nn (k;P)}, (9) 

and f nn , its leptonic decay constant 

U n & P, = Z 2 tr f A ir* 75 7^ S(q + )Vi n (q- P)S(qJ) , (10) 

Jq 

where the trace is over colour, flavour and spinor indices. The lowest mass pole con- 
tribution, denoted by n = 0, is the ground state pion, which receives the spectroscopic 
assignment N 2s+1 Lj = 1 in the naive qq quark model. The n > 1 pseudoscalar meson 
poles correspond to those excited states of the pion that would receive the spectroscopic 
assignments (n + 1) in the quark model. 

The pseudoscalar vertex, Tl(k; P), which appears in Eq. (J3J), also exhibits such a pole: 

^ P) U mL ,0 = termS + P 2 P + m 2 T iS k -> P ) > ( H ) 

ip Wn (0 5 ij = Z,tr [ A ^ l5 S(q + )Ti n (q; P)S(qJ) . (12) 

Jq 

In QCD it therefore follows as a necessary consequence of chiral symmetry and its 
dynamical breaking that for any pseudoscalar u-d meson [ ITT] 

U n ml n = [m u {Q+m d {Q]p« n {Q- (13) 

(The generalisation to SU(Nf)-fi&vom is straightforward, following Ref. [HE].) The so- 
called Gell-Mann-Oakes-Renner relation for the ground state pion appears as a corollary 
of Eq. (fT3~]) [HHj, and another important corollary of Eq. (|T5|l. valid for pseudoscalar mesons 
containing at least one heavy-quark, is described in Ref. [EOJ. In addition, the proof of 
Eq. ()13j) establishes [El that in the chiral limit 

/° v := lim U M = Q; (14) 
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i.e., in the chiral limit the leptonic decay constant vanishes for every pseudoscalar meson 
except the pion. This exact result is a constraint on models and methods used to search 
for exotics and hybrids. 

The rainbow-ladder truncation obtained with the interaction of Eq. (j7J) provides an 
excellent description of ground state pseudoscalar mesons [ |2] . Nevertheless, the long- 
range part of this interaction expresses a model and its veracity can be tested further by 
studying excited pseudoscalar mesons. We have begun this process, focusing on the n = 1 
pseudoscalar meson (2 in the quark model). 

In our studies to date we have only retained the i^E ni (k; P) piece in Eq. (jUJ). Employing 
this expedient for the ground state pion leads to a 25% underestimate of ra^ . Moreover, 
in neglecting the pseudovector components, we omit those terms which in the pion's 
rest frame signal the presence of quark orbital angular momentum. In QCD they are 
necessarily nonzero [ [EH E] . 

To determine m 7rn=1 we first solved a rainbow truncation of the gap equation; viz., 

Sip)' 1 = Z 2 (i 1 - V + m bare ) + jf * 47ra((p - q) 2 ) D^(p - q)^ lfl S(q)^ lu , (15) 

where D { ™ e (k) is the undressed gauge boson propagator. This allows one to complete the 
specification of the ladder-truncation Bethe-Salpeter equation's kernel: 

r„ n (k; P) + J\na((k - qf) D^ c (k - <?)y 7 „<%+) T^q; P) S(q_) y 7 „ = , (16) 

which can be solved as described in Ref. [Ej to obtain the mass of the ground state and 
excited pseudoscalar mesons. 

One should note that in relativistic quantum field theory the scalar functions in Eq. (J§|) 
depend on two variables: k 2 , k ■ P. In quantum mechanics k ■ P = 0, of course, because 
of "on-shell" constraints on the constituents, and this degree of freedom is absent. As a 
matter of practice, we solve Eq. (|16j) by using a Chebyshev decomposition; i.e., we write 2 

■*V max 

E nn (k 2 ,k-P;P 2 )= ^„(£; 2 ;P 2 m(cos/3), (17) 

i=0,2,4,... 

where {Ui(x); i = 0, . . . , oo} are Chebyshev polynomials of the second kind, and increase 
iV max until m 7Tn stabilises. This procedure has the merit of limiting the amount of computer 
memory required. It also converges very rapidly for the ground state pion, with only the 
leading Chebyshev moment being required [QB], as apparent in Fig.|3J 

We have obtained the mass and amplitude for the ground state pion, using the complete 
expression in Eq. ©, and found: f no = 0.092 GeV; m wo = 0.14 GeV; p wo = (0.81 GeV) 2 , 
at a current-quark mass m<j(l GeV) = m u (l GeV) = 5.5 MeV, reproducing the results in 
Ref. [E]. We can report in addition that using solely the i^ 5 E m {k\ P) term in Eq. (jUJ) 

m ni = 1.1 GeV , (18) 

cf. m^ pt ' = 1.3±0.1 GeV. In omitting the pseudoscalar meson's pseudovector components 
one violates Eq. (0J) [ E] and hence we are not yet in a position to verify that Eq. (|14|) is 

2 The Bethe-Salpeter amplitude for charge conjugation eigenstates is even under fc • P — > —k ■ P and hence 
only even Chebyshev moments contribute. 
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Figure 3. Chebyshev moments, Eq. (JT7J), of E^ {k;P) = E gr (k; P) and E ni {k;P) = 
E ex (k;P). The leading moment is sufficient to accurately represent the ground state 
while a veracious description of the first excited state requires inclusion of the second mo- 
ment. Note the appearance of a single node in each of the first excited state's Chebyshev 
moments. (The amplitudes are canonically normalised.) 



preserved in the rainbow-ladder truncation. Nevertheless, we estimate f Wn=1 < 1.5 MeV 
at the physical current-quark mass [ITTj. 

We plot the leading Chebyshev moments of E ni (k; P) in Fig.0 wherefrom it appears 
that the first two moments are sufficient for an accurate representation of this amplitude 
but confirmation must await the inclusion of the pseudovector components. It is worth 
remarking that the moments of E ni (k; P) each possess a single node and exhibit tails that 
extend to larger values of k 2 than their analogue in the ground state. These observations 
are redolent of those one could make in comparing Fourier transforms of the radial wave 
functions for bound states in a potential well. This fact emphasises the intuitive character 
of Bethe-Salpeter amplitudes. 



4. EPILOGUE 

There are many additional applications of interest to this community, among them: the 
ab initio calculation of electroweak and transition form factors for other mesons, and titt 
scattering [122]; and a calculation of the pion's valence-quark distribution function, whose 
discrepancy with extant data raises difficult questions [ 23j. A pressing contemporary 
challenge is the extension of the framework to the calculation of baryon observables, 
aspects of which are beginning to be understood [123]. 
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